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Abstract: We show that the two couple equations derived by approximate symmetry 
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1 Introduction 

Group theoretical methods based on local and nonlocal symmetries provide remarkable tech- 
niques to analyze partial differential equations (PDEs) such as constructing group invariant 
solutions, studying integrability etc., and also build a ladder to explore profound understand- 
ings of the physics of the underlying problems such as detecting conservation laws [TJ[2] . The 
prerequisite of achieving these goals is that the PDEs under study possess affluent symme- 
tries, nevertheless, for the perturbed PDEs (containing a small perturbed parameter), the 
classical symmetry method losses its superiorities because the small parameter disturbs the 
symmetry group properties of the unperturbed equation, even makes the studied PDEs hold 
few symmetries. 

In view of the drawbacks of classical symmetry method for conducting this type of equa- 
tions, the integration of perturbation analysis and symmetry group theory makes two well- 
known approximate symmetry methods (ASM) emerge [31 Q]. The originators of the two 
approaches are Baikov et.al. [3j and Fushchich and Shtelen [4] who employ the standard 
perturbation techniques on the infinitesimal operators and dependent variables respectively. 
The comparisons of the two ASM which are illustrated by several specific equations may 
be found in [5j[6]. Quite recently, the extended combination of homotopy conception with 
perturbation technique and symmetry method produces approximate homotopy symmetry 
method (AHSM), which is suitable for studying the problems which don't contain small 
parameters (7j- 

The distinct merit of AHSM is that the accuracy of approximate homotopy series so- 
lutions can be controlled by adjusting the convergence-control parameter, however, it is at 
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the expense of the greater computation than ASM [8]. On the other hand, we find that 
the frameworks of ASM and AHSM for the perturbed PDEs are quite similar except for the 
different studied objects (see subsection 2.1 below). Hence, a natural question occurs: what 
is the relationship between ASM and AHSM, does there exist a transformation or others to 
link them? In [8] , the authors study this problem which is exemplified by a class of perturbed 
nonlinear wave equations and show that first-order coupled equations and higher-order simi- 
larity reduced equations are connected by two scaling transformations respectively. However, 
the general results for the connections between ASM and AHSM for general perturbed PDEs 
are still open. 

The purpose of our paper is to give a positive answer about this question. We show that 
the two couple equations obtained by ASM and AHSM are equivalent under a transformation. 
Based on this transformation, one can construct approximate homotopy series solutions 
directly from the known approximate solutions by ASM without performing AHSM for the 
governing equations. 

The outline of the paper is as follows: In section 2, we briefly review some basic notions 
and give the main results. Section 3 concentrate on the applications of the results to the 
Cahn-Hilliard equation. The last section contains conclusion and discussion of our results. 

2 Basic notions and principles 

In this section, we first briefly review the main ideas of ASM and AHSM, and then give the 
results about the relationship between the two approaches. 

2.1 Basic notions 

We take the following /cth-order perturbed nonlinear PDE 

E = E (x, u, u {1) u {k) ) + e Ex(x, u, , . . . , = (1) 

into account to illustrate the main ideas of ASM and AHSM, throughout the section, with 
independent variables x = {x\, . . . , x n )(n > 2), a single dependent variable u, vfj> = 

d j u/dx 3 ^ . . . dx% with \-j n = j (j = 1, . . . , k) is a collection of all jth-order 

partial derivatives, e is a small perturbed parameter. Note that we only use one variable 
u to keep the notation simple, the generalization to systems being fairly obvious. Also 
Ei = Ei(x, u, . . . , u^>){i = 0, 1) if no special notes are added. 

Let's first recall the definition of classical symmetry. Consider the infinitesimal generator 

n 

x = Y1 n °> Ul ■ ■ ■ ' u{k) ) d ^ + 7 ?( x ' u o,ui..., u {k) )d u , (2) 

where and r\ are smooth functions of their arguments. Let be the fcth-prolongation 
of the generator X calculated by the well-known prolongation formulae [U[2]. 

Definition 1. (Symmetry or exact symmetry [HE]) The Lie symmetry X of the form 
([2]) is admitted by Eq. (TTJ if and only if 

X^E\ E=0 = 0, 
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where the notation \e=o means evaluated on the solution manifold of E = 0. 

The ASM originated by Fushchich and Shtelen employs a perturbation of dependent vari- 
able and then approximate symmetry of original equation is defined to the exact symmetry 
of the system corresponding to each order in the small parameter. 

Specifically, expanding the dependent variable with respect to the small parameter e 
yields 



u 



1=0 



where ui = ui(x) (I = 0, 1, 2, ... ) stand for unknown functions throughout the paper, then 
substituting it into Eq.([TJ and vanishing the coefficients of all different powers of e, one 
obtains a coupled system 



Eo\ e= o — 0, 
+ i 



a 



-i-l 



|e=0 



0. 



1,2, 



(3) 



Definition 2. (Approximate symmetry [4]) The iih order approximate symmetry of Eq.(pQ) 
is defined to the exact symmetry of the first % + 1 equations in Eq.([3|). 

For approximate homotopy symmetry of Eq.(pQ), one generally consider the following 
simple homotopy model 



H(u,i 



I - 0q)E o + qe(l - 0)E l = 0, 



(4) 



where, hereinafter, q G [0, 1] is an embedding homotopy parameter and 9 denotes the 
convergence-control parameter. The above homotopy model has the property H(u, 0) = 
E , H{u, 1) = E. 

Assuming Eq.(pQ) has the homotopy series solutions in the form 



(5) 



1=0 



where u\ solve 



Eo\ q= o — 0, 



q 



dqi 



Oq 



i-l 



+ ie(l 



d l - l E l 



Oq 



i-l 



1 9=0 



0. 



1,2, 



(6) 



Equivalently, we rearrange Eq.([6]) by eliminating d' 1 Eo/dq t 1 in the (i + l)th equation 
by means of the first i equations and transform it into the following form 



E, 



0\q=0 



0. 



- dq l 



i-l 



ii—1—k 



k=0 



i\ d k Ef 
k\ dq k . 



1 9=0 



1,2,.... 



(7) 
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Definition 3. (Approximate homotopy symmetry [7]) The ith order approximate ho- 
motopy symmetry of Eq.([T]) associated with homotopy model ((H) corresponds to the exact 
symmetry of the first i + 1 equations in Eq. (J6|) or Eq. (J7|) . 

Obviously, the first equation of Eq.Q is the same as the first one of Eq.©. When 9 = 0, 
Eq.© is equivalent to Eq.(J7]) under the scaling transformation Ui = e l Ui (i = 0,1,2,...) 
owing to {d' L Eq/ dq % )\ q=Q is linear about Ui and its derivatives (see Lemma 1 below). 
However, the parameter 9 is used to adjust the convergence of homotopy series solution 
and may not be zero, thus it is the intention below to search for the connections under the 
condition 9^0. 



2.2 Main results 

In what follows, we first prove two basic lemmas and then establish the connections between 
Eq.® and Eq.©. 

Lemma 1. (d n Eo/dq n )\ q ^ is linear about u n and its derivatives \ where n(> 1) are 
positive integers and A = 0, 1, . . . , k. 

Proof: We take u^°\= u),vP~\ . . . of Eq as independent arguments and expand u 
as ([5]), then by means of the generalized Faa di bruno's formula [11], nth-order derivative of 
Eq with respect to q is 



d n E \ \- nl d r E 

8q n _ o 1 ' n TT C1\r TT W I X <9[u(°) ]Pl C?^ 1 ) . . . d[u^] Pk + 1 



i=l i=l j=l 



X 

1=1 



where the respective sums are over all nonnegative integer solutions of the Diophantine 
equations, as follows 

■ ri + 2r2 H \-nr n = n, 

o 

^ ■ Pml + Pm2 H h Pm(fc+1) = r m, m = 1, 2 



, ... ,71, 



«* = pi + P2 H h Pfc+i = ri + r 2 H h r n , 

Pj=Plj+P2jH hPnj, j = 1,2, . . . , k + 1, (9) 



and 



^E(tV-^>, ^0,!,...,, (1 0) 
Obviously, when (/ = 0, all terms in Eq. (|10p vanish except for u\ X \ Similarly, Eq.([8|) with 
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q = becomes 
d n E 
dq n 



EE-E- 



1 9=0 1 



n 

*n[» 



fc+i 

n n n 

i=i j=i 



^ 

a^w^a^w^ . . . 9[«(fc)]Pfc+i 



1 9=0 



(0) 



i=l 



Pil 



(1) 



Pi 2 



(k) 



Pi{k+1) 



(11) 



In Eq. tjlip . only when i = n, the terms containing u n ,(or its derivatives) occur. At the 
moment, the solution of the first equation in Eq.Q r n = 1, rj = 0(i 7^ n) generates the terms 
uniquely, then other equations give r = 1 and p» = = 0(s 7^ i) and p n j = l,p n s = 0(s 7^ 
i). Hence, the arguments u n and its derivatives only occur in the following terms in Eq.ljlip 



dE, 



dE 



,(1) 



d£ 



u 



(fc) 



which are linear about u n and its derivatives since (8Eq/ du^)\ _ (i = 0, 1, . . . , fc) are only 
functions of x, iio and their derivatives. Meanwhile, all other terms with q = are only func- 
tions of x, u, U\, . . . , u n -\ and their derivatives, thus Lemma 1 follows. The proof ends. □ 
Lemma 2. Substituting the transformation 

i-l 



'0 



' 



(A) 



E 

i=o 



^ ^V^-, « = 1,2,..., A = 0,1,. 



, k. 



into (d n Eo/dq n )\ q=0 , after deleting the symbol^, [3 n i?o(0 , u^)/dq n ]\ is equal to 



d n E { 



- dq 



^0 y 

n / j 



i=2 



(i - 1)! - 2 



n-l\ „ n _ i+1 &^Eo 



dq 



i-l 



1 9=0 



(13) 



where [d n Eo(6 : u, . . . , vS ') / dq n ]\ q=Q denotes the transformed expression via (| 12j) and 
{d l Eo/dq l )\ q=0 (i = 2, 3, . . . , n) stand for the untransformed ones. 

Proof: The main idea of the proof is to regard [d n Eo(0 : u, . . . , vS k ')/dq n \\ =0 as a poly- 
nomial of 9 and search for the coefficients of different degree of 6 in it. 

Specifically, transformation (|12p makes Eq. (jllj) become 



d n E 



EE-E 



nl 



1 



n 

i=l 



n k+1 

n n p«i 

i=ij=i 



9=0 



E 



i=o 



j-'ii 



i-l 

E 



EE'E 



n 



fc+i 

' = 1.7 = 1 



eJu i-j 



d r E 



Pi(k+1) 



19=0 



x[2<°) + ...+ 



+ 



p-l^( fe )lp„(fe+l) 



(14) 
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In particular, 



j=0 



u 



(A) 
i-j 



E 



a v0 ! • • • 



E 



a v0 \ . . . o„(i_i)! 

ia J/ i+2a„2H h(i— l)<V(i-i) 



y(i-i) 



(A) 



(15) 



where A = 0, 1, . . . , k and the sum is over all nonnegative integer solutions of the Diophantine 

equation a u0 + «Vi H 1- = Viv 

For Eq. ()14p . we first consider the terms without containing parameter 0. That is to say, 
let = in Eq.flHJ), we nave 



d n E n 



Q a n L^i L^i n k+1 

" nn py\ 



0[u(°)]pi0[u(i)]P2 ...0[&W]pm 



1 9=0 



i=l J=l 



7T, 

«n[» 

i=l 



(0) 



(1) 



P( 2 



u 



(k) 



Pi{k+1) 



(16) 



which has the same form as Eq.(|lip. 

Secondly, we look for the terms only involving in Eq. (|14p . The general formula of this 
type of term appears as follows 



2,3,... 



3=0 



J[^]PiU+D x ■ ■ ■ x f j [u^fw- 1 ^^ x • • • x [^')]P»y+D, m 



At this time, the parameters satisfy 

^ : ri H h (m - l)(r m _i + 1) + m(r m — 1) H \-nr n = n — l, 

o 

X] : P(m-l)l + P(m-1)2 H 1" P(m-l)(fc+l) = r m-l + 1, m = 2, . . . , n, 

m— 1 

^ : Prnl + Pm2 H h Pm(*+1) = r m - 1, m = 2, . . . , n, 

^ : + p S 2 H h p s (fc+i) =r s , s = 1, . . . , m - 2, m + 1, . . . , n, 

s 

r = pi + p 2 H h Pfc+i = r\ + r 2 H h r„, 

Pi = Pij + P2j H h Pnj, j = 1, 2, . . . , k + 1. 



(17) 



Especially, due to r±, . . . , r n are nonnegative integers, thus from the first equation of (H 



(3 



one has r n = 0, which make condition (|17j) become 

■ n-\ h (m - l)(r m _i + 1) + m(r m — 1) H h (n - 1) r n _i = n - 1, 

o 

X] : + P(m-1)2 H 1" P(m-l)(fc+l) = r m-l + ^ 771 = 2, . . . , n - 1, 

m— 1 

X] : Pml + Pm2 H h Pm(fe+1) = r m ~ 1, 771 = 2, . . . , 77 - 1, 

m 

p s i+p S 2H ^P s (fc+i) = r s , s = l,...,m- 2,m + l,..., n, 

s 

r = pi + p 2 H h Pfc+l = 7*1 + 7*2 H h 7* n -l, 

Pj =Plj +P2j H l"P(n-l)j, i = 1,2, + (18) 

The condition f)18|) is just the requirement of the parameters in (d n Eo/dq n )\ _ , and 
the coefficient of 9 in Eq. (fT3|) corresponds to (d n Eq/ dq n )\ q=0 multiplied by 77(77 — 1). 
Hence, we prove that Eq.flHJ) involves 77(77 — \){d n ~ l Eq / dq n ~ 1 )\ =0 . 



Enlarging the degree of 9 and repeating similar procedure as above up to 



qn-2 



WC 



find that Eq. (|14p contains 



rfl-l 



i=2 



\q=0 



, thus at last, we consider the 



maximal degree of 9, that is 9 n ~ 1 , which only occurs in the terms containing u n or its 
derivatives. Hence, this case corresponds to r = 1 in Lemma 1, which implies the terms 
corresponding to 9 n ~ 1 in Eq.flTlJ) reads 



77! 



dE (o) dE (i) dE 

77] + - — U) + • • • + - ; 77{ 



9u( fc )| 9 =o 1 



(19) 



which is just the (dEo/dq)\ multiplied by 77 !. This prove it. □ 

Next, for 9 7^ 0, we give the main theorem about the connections between Eq.([3D and 
Eq.©. 

Theorem 1. Eq.([7|) is equivalent to Eq.([3|) under the transformation 



1-1 



I - 1 



^[ e (l-0)]^g, 



(20) 



where Z = 1, 2, . . . , A = 0, 1, . . . , k. 

Proof. We start with Eq. ([7]) and list it for different order of q up to I in detail as follows 

Eo\ q= o = 0, 

dE 



- dq 

d 2 E 



+ e(l-9)E 1 



\q=0 

fdEx 



+ 2e(l-0) — - + 0E 



dq 2 ' V dq 



\q=0 



0. 



oq z oq 



\q=0 



0. 



(21) 



^+Ki-«)(^ + (i-D^ + - + a-D"^ 



9=0 
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The proof of Lemma 1 tells us that d l Eo/dq l with q = is linear about u\ (I = 1, 2, . . . ) 
and its derivatives, thus we first consider transformation (|12p for Eq. (|21l) . 

The first equation in Eq, (|2ip is unperturbed equation which is the same as the first 
equation in Eq.Q. For the second equation, by Lemma 1, (dEo/dq)\ =0 is linear about U\ 
and its derivatives, thus one can adopt a scaling transformation to convert it to the same 
form as the second one in Eq.©. 

Next, we consider the third equation in Eq. (|2ip . Take the first three transformations 
of (|12p into the third equation in Eq. (|2ip and minus the second equation multiplied by 26, 
then by means of Lemma 2, one obtains 



rj2 Q 

g^E (x, u . . . , + 2e(l - 0)—^ (x,u..., u^] 



which has the same form as the third equation in Eq.([3]) 



(22) 



1 9=0 



d 2 E 2 dE x 



Oc 



Oc 



\e=0 







except for the coefficient of dE\/dq. 

In the same way, substituting the first four transformations of (|12p into the fourth equa- 
tion and minus the sum of the second equation multiplied by 6 6 2 and Eq. (|22p multiplied by 
6$, also using Lemma 2, we have 



r 8 3 8 2 
q^E (x, u . . . , v^) + 3e(l - 0)-^E l (x,u..., u^] 



which has the same form as the fourth equation in Eq.([3|) 

^2 ; 



0, 



(23) 



1 9=0 



d 3 E ^E, 

~T~ O- 



de 3 



Oc 







except for the coefficient of d 2 E\jdq 2 . 

Enlarging / inductively and assuming that the first I equations of Eq. ([2ip have been 
transformed to 







t-i 



^E (x,u...M k) )+e(l-l)(l 



.dq l 



& 



1-2 



dq l 



; Ei(x,u. 



0, (24) 



9=0 



then insert the first (I + 1) transformations of (|12p into the (I + l)th equation and minus the 
sum of the ith equation of the form ([24]) multiplied by (jzjyr ('l^)^'^ 1 (* = 2, 3, ... ,/), we 
convert the last equation of Eq. (|2ip to the form 



^jE (x,u..., ) + d(l- 0)-^ T E 1 (x,u..., ) 

which has the same form as the (I + l)th equation in Eq.© 

-d l E n , d l ~ l E 



(25) 



9 = 



+ 1 



Q € l-1 
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except for the coefficient of d 1 E±/dq . 

At last, by Lemma 1, using the scaling transformation 

4 A) = 4 A) , u[ X) = e(l - 9)u[ X \ u {X) = e(l - 0)4 A) , . . . , uf A) = [e(l - 6)] l u {X \ (26) 

one can convert Eq. (|25p to the (/ + l)th equation of Eq.© by ASM. 

In summary, combining transformations (112j) and (126p . we obtain transformation (120p . 
which connects Eq.([3]) and Eq.([7D. This proves it. □ 

Transformation ()20|) in Theorem 1 build a bridge between AHSM and ASM for per- 
turbed PDEs, that is to say, there is no need to perform AHSM to the governing equations 
because the desired approximate homotopy series solutions can be obtained by acting the 
transformation (|20p on the approximate solutions by ASM. 



3 Applications 

In this section, we apply transformation (|20p to study the connections between ASM and 
AHSM for the Cahn-Hilliard equation piTO] 



ut = - [F(u)u x ] x - eu xxxx , (27) 

which is one of the simplest models to characterize the process in the context of the continuum 
theory of phase transitions. In Eq J27p . x, t are two independent variables and u is a scalar 
dependent one, the nonlinearity F(u) is the derivative of a double-well potential with wells 
of equal depth and < e < 1 measures the thickness of an interface separating the two 
preferred states of the system. 

3.1 The connections between the couple equations 

From Eq. (|27p . one has 

E (u) =u t + [F(u)u x ] x , E\{u) = u xxxx . 

By means of the general formulas of couple system ([7|) , we get the following couple equations 
by AHSM 



,t + [ E F r] {U ° ] -A < u h ■ ■ ■ 4>k-i, x ] + e(l ~0)J2 ek ~ 1 ~ i ^ xxx = 0. (28) 

i=0 "' '"' 1 x i=0 

where, hereinafter, n_i = 0, uj* = (u lt ) ji , u k>t = du k /dt, u itXXXX = d A Ui/dx 4 , j +iiH H 

ji = J) ^oio + hji + • • • + hji = i, < j < i and Zo, h, . . . , 1% are not equal to zero and mutual 
inequivalent. k,ji,j,k (i,j,k = 0,1,2,...) are nonnegative integers and satisfy the above 
relations. 

Now, we use transformation (|20p to convert Eq. (j28j) to the couple equations by ASM. In 
order to give specific procedure about how transformation (|20p works for Eq. (|28p . we first 



k-l 
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multiply both sides of Eq. (|28p with k\ in order to make it with integer coefficients and then 
list it for different order of q in detail as follows 

uo,t + [F(uo)u 0>x ] x = 0, 

u±,t + [F(u )ui\ X x + e(l - 9)u , xxxx = 0, 

2u %t + [F'{u )u\ + 2F(u )u 2 } xx + 2e(l - 9)[ U\ xxxx + ^Uo xira] — 0, 
6n 3ii + [F"{u Q )u\ + QF'{uq) Ui u 2 + 6F(u )u 3 ] xx 

,xxxx 

+ 0ui 

,xxxx ,xxxx] = 0, 

(29) 

U k\F^(uo) 



-l-i. 



kl «*.« + [ E fSx^X ■ ■ ■ , + < l - °) kl £ Qk 

Like the proof of Theorem 1, we first consider transformation 

S o, 4 A) = £ ~ ^^i-i- fc = l,2,..., A = 0,l,...,fc. (30) 
i=o V z ' 



»0 



By means of Theorem 1, the first two equations of Eq. (|29j) are equivalent to the cor- 
responding ones by ASM under a scaling transformation uo = uq,u\ = e(l — 9)u±, thus 
we consider the third equation. Take the first three transformations of (130p into the third 
equation and minus the second equation multiplied by 29, one obtains 

2u 2 ,t + [F'{uq)u\ + 2u 2 F(u )] xx + 2e(l - 9)u 1)XXXX = 0. (31) 

Similarly, substitute the first four transformations of (|30p into the fourth equation and 
minus the sum of the second equation multiplied by 69 2 and Eq. (j3~Tj) multiplied by 69, we 
have 

6u 3)t + [F"(u )ul + 6u 2 u 1 F'(uq) + 6u 3 F(u )] xx + 6e(l - 9)u 2jXXXX = 0. 

Suppose the first k equations in Eq. (f29j) be in the form 

fc-i 



( k-l)\F^(u ) j0 n k 

- - ■ U,U, . . . UjUk-l-i a 

j \j!\...ji\ h l * 



(fc-l)!« fc _i )t + [E 

+e(l-9)(k-l)\u k _ 2tXXXX = 0, (32) 



then insert the first (k + 1) transformations into the (k+ l)th equation and minus the sum of 

kl - 



the zth equation of the form (f32j) multiplied by (^_2)^ fc * +1 = 2, 3, . . . , fc), we convert 



the last equation of Eq. ([29j) to the form 
k\F®( 

i=0 



k\ U k t + [E ^ U °\ u ^ U h ■ ■ ■ U l u k-i.x] + e(l - Uk-l,xxxx = 0. 



Finally, since the ith equation of Eq. (|29p is linear about itj and its derivatives (i 
1,2, ... , fc), thus integrating the scaling transformation 

u Q —u , u x — 6(1 — V)u x , u 2 —e(L-V)u 2 , ...,u k —[e[L — U)\u k , 
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and (I30p . one obtains (|20l) . which convert Eq. (l28p to the form 



:0 J0!jl! • • • Jj! 



■f- /c. Uk—l,xxxx — 0, 



after deleting the symbol , which is just the couple equations via formulas ([3]) by ASM. 
3.2 Homotopy series solutions through the transformation 



In this subsection, we use transformation (|20p to obtain approximate homotopy series solu- 
tions from the approximate solutionsby ASM. 
I. Traveling wave form 

By ASM, the traveling wave solution with third-order precision for F(u) = 1/uis given 
by u = uq + eu\ + e 2 U2 + e 3 us after selecting proper parameters, where 

1 3 774 51273 

""0 = -?— : T> u l = 7T, 71' u l = — o , 77i u 3 



a(at-x)'^' a 2 (at-xf 10a 3 (at -x) 7 ' 10a 4 (at - x) 10 ' 

and a is wave speed. 

By virtue of the first four transformations of (|20p . i.e., 

Uq = uo,«i = e(l - 9)u\, 

u 2 = e(l-9)[e(l-O)u 2 + 0u 1 ], 

u 3 = e(l - 9)[e 2 (l - 9) 2 u 3 + 20e(l - 0)5 2 + # 2 5i], 

we obtain a third-order approximate homotopy series solution of Eq. (j27p in the form u 
uq + qui + g 2 U2 + q 3 U3, which is 

1 3eq(9 - 1) 3eg 2 (# - l)[5a0(at - x) 3 + 129e(# - 1)] 

U — — 7T. I iT7 777 "4" 



a 2 t — ax a 2 (x — at) A 5a 3 (at — x) 7 

3eq 3 (6 - l)[10a 2 fl 2 (x - at) 6 + 5We9a(9 - l)(at - x) 3 + 17091e 2 (# - l) 2 ] 

10a 4 (x - at) 



' 1 TTZITZ ITuo ' '33) 



Inserting (j33j) into Eq. (|27p and putting a; = l,t = 0.1, e = 0.01, we get figure 1 which 
demonstrates the error of solution (]33p with respect to the auxiliary parameter 9. 

From Figure 1, in particular, assuming convergence-control parameter 9 = 0.5478, one 

has 

u t + [F(u)u x ] x + eu xxxx « 4.70 x 10~ 6 . 

II. F = u 

We construct the following approximate solutions with fourth-order precision after as- 
signing proper value to some parameters by ASM 

x 2 e 30e 2 t 46440e 3 i 2 804646440e 4 t 3 



u = — + —s 



U x 2 x 6 7x 10 203x 14 
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0.1 




Figure 1: x = 1, t = 0.1, e = 0.01, horizontal axis denote 8, vertical axis denote the error 



By virtue of the first five transformations of (|20p . a fourth-order approximate homotopy 
series solution follows 



x 
+ 



eq(9-l) eq 2 (9 -l)(390e(9 -l)t + 139z 



x 2 13x 6 
eq 3 (9 - 1) (603720e 2 (6» - l)H 2 - 5460e6»(6» - l)tx 4 - 9W 2 x s ) 



+eq\6 - 1) 



91x 10 

-10460403720e 3 ((9 - l) 3 t 3 + 52523640e 2 6(9 - l) 2 t 2 x 4 



2639x 



ii 



237510e6» 2 ((9 - l)tx 8 + 26399 3 x 12 



(34) 

2639x 14 J v ; 

When the constants and variables in solution (|34|) are specified byx = 1, t = 0.1, e = 0.01, 
the relationship between the error of fourth-order precision approximate homotopy series 
solutions of (|34|) and parameter 9 is sketched in Figure 2. 




Figure 2: :r = l,f = 0.1,e = 0.01, horizontal axis denote 6, vertical axis denote the error 

Especially, taking convergence-control parameter 9 = 0.7015, one has 

u t + [F(u)u x ] x + eu xxxx w 1.59 x 10~ 6 . 
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4 Conclusion and discussion 



In this paper, we show that the two couple equations derived by ASM and AHSM are con- 
nected by the transformation (|20|) . As a result, approximate homotopy series solutions can 
be obtained by acting the transformation on the known solutions by ASM. As an applica- 
tion, the Cahn-Hilliard equation is considered to examine the results. This transformation 
may exert an important role in studying the perturbed PDEs because it makes researchers 
facilitate to study the results by AHSM through the ones by ASM. 

In addition, it is worthwhile pointing out that, for Eq. (|27p . the unique pair of different 
operators Y\,X\, i.e., the approximate homotopy symmetry operator 

Yi = --xd x + \- iUi l W u i~i]d Ul 
i=i 

and the corresponding approximate symmetry one 

Xi = --xd x + iuidzt 

are equivalent under the transformation 

t=i (I - 1\ I - i M „ 
u = u , u i = Z^\ ■ \—j- eu i-h 1 = 1,2,..., 

other than (j'20|) . Hence, the relations between the symmetries of two equivalent equations 
may form an interesting topic. 
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